Spin and charge excitations in a Three-Legs Fermionic Ladder: a 
Renormalization-Group study. 
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We study the spin and charge phase diagram of a three-legs ladder (at zero temperature) as a 
function of fermion density and of transverse single-particle hopping by means of a Renormalization- 
Group analysis rigorously controlled in the weak-coupling limit. 

Periodic boundary conditions in the direction transverse to the ladder produce frustrated mag- 
netic excitations yielding a spin-gapped phase in a large region about and at half filling. Spin 
correlations are instead enhanced when open transverse boundary conditions are considered, yield- 
ing an ungapped phase in a wide region about half filling (as observed in the Sr n _iCu n +i02n stripes 
compounds) and at half filling down to a critical value of the transverse hopping t j_ . At that critical 
value, the system undergoes a Mott transition at half filling by decreasing t±. 



PACS numbers : 72.15.Nj, 64.60.Ak, 71.27. +a 

The problem of coupled one-dimensional (1-D) systems 
of interacting fermions (ladders) |j]-|4| has recently re- 
ceived considerable attention among condensed-matter 
physicists due to the possible connections with high-Tc 
superconductivity. In addition, these systems present un- 
usual properties whose analysis can give valuable new 
insights into the properties of low-dimensional quantum 
systems. An interesting topic addressed in this field is 
the contrast between the magnetic properties of ladders 
with even and odd numbers of legs, in analogy with 1-D 
systems with integer and half-integer spin. |^,|5| 

This theoretical interest has been further prompted by 
the discovery of the Sr„_iCu n +i02 n compounds, whose 
lattice structure shows strips of (n + l)/2-legs fermionic 
ladders. |(| Different measures on these systems, show 
that a spin gap (signaled by an exponential falloff of the 
magnetic susceptibility as the temperature decreases) oc- 
curs in the even-legs compounds (at least at half filling), 
while the odd-legs systems remain ungapped. 

In a strictly 1-D spin 1/2 system the low-lying spin and 
charge excitations are completely separated and the asso- 
ciated correlation functions show long-range oscillations 
with an anomalous power-law decay, whenever high- 
momentum-transferring processes (backward and umk- 
lapp scattering) can be neglected, allowing for the low- 
energy problem to be solved exactly. R] At half filling, 
however, the coupling associated with the umklapp pro- 
cess becomes relevant and opens a gap in the charge exci- 
tations causing the system to be an insulator (at T = 0) 
for arbitrary small repulsion. Q This does not happen for 
the repulsive backward scattering which remains irrele- 
vant leading to quasi-long-range-ordered spin degrees of 
freedom. The problem is thus to understand when a given 
processes is important (relevant) for the low-energy (or, 
equivalently, long-range) properties of the system. This 
issue is best addressed in the framework of the renormal- 
ization group (RG), whereby one follows the evolution 
(flow) of the couplings as the system is observed by larger 
and larger length scales. ||] 

The same question can be addressed for multi-legs sys- 
tems, where the analysis is complicated by the larger 
number of parameters. Extensive RG studies have been 
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carried out for two-legs ladders, [1 
terns with an higher number of legs 
these calculations clearly show that the crossover from a 
1-D system to a system with many legs is far from be- 
ing smooth at least concerning spin excitations, which 
are exponentially suppressed in even-legs system, while 
they are quasi-long-range ordered in systems with an odd 
number of legs. Data from the Sr ra _ 1 Cu„ + i0 2rl ladder 
systems agree with these results. [0] A smooth behav- 
ior as a function of the number of legs should be then 
recovered when either even-legs or odd-legs systems are 
considered separately. In this sense, a three-legs system 
is the first continuous extension of the purely 1-D system. 
Anomalous exponents for the uniform charge degrees of 
freedom should instead continuously decrease with the 
number of legs, at least away from half filling. [Q 

It is thus worthwhile to address the properties of an 
odd-legs system, by extending the analysis also to the 
charge degrees of freedom and by considering the behav- 
ior of the system upon doping. In this paper, we study 
the phase diagram of a three-legs spin 1/2 fermionic lad- 
der by means of weak-coupling RG approach like the one 
applied to one || and two coupled one-dimensional 
systems, here restricting to one-loop order. We mainly 
analyze the case of open boundary conditions in the y di- 
rection (OBC), since it represents the experimental situ- 
ation and it is free from frustration in the magnetic inter- 
actions. Our results yield indeed ungapped spin modes in 
a wide region about half filling, and at half filling for an 
almost isotropic system. In the same region, the system 
is a conductor, whereas a metal-to-insulator transition 
appears at half filling for low values of t p . 

Specifically, we consider a system of 3 identical one- 
dimensional single-band models (legs), each with band 
energy e(k) (\k\ < n), coupled via a single-particle trans- 
verse hopping t p and supplemented with a generic short- 
range interaction. For the sake of definiteness, we will 
specialize to the Hubbard model with on-site repulsion U 
and nearest-neighbor hopping t, although we will even- 
tually allow for extensions. By diagonalizing the single- 
particle term containing t p and e(k) one obtains 3 "par- 
allel" bands with energies e(fc, a) = e(k) — T p (a) where 



a = 1, . . . , 3 labels the 3 bands, and T p {a) — \f2t p (a — 2) 
or 2 t p cos 2,r (^~ 1 ) for OBC and periodic boundary con- 
ditions in the y direction (PBC), respectively. 

Within an RG approach one first "integrates out" the 
degrees of freedom far from the Fermi surface (consisting 
of 6 points when all bands cross the Fermi surface) and 
restricts the range of the momenta to regions up to a 
cutoff A about each Fermi point, such that, within each 
of these disconnected regions, the dispersion can be ap- 
proximated with a linear function and the interactions 
can be taken as fc-independent. The renormalization of 
the parameters (Fermi momenta, Fermi velocities, and 
interactions) during this "preliminary" integration can 
be neglected in the weak-coupling limit and their bare 
values can be used as "input" for the linearized system. 
I 

The generic interaction limited to the restricted regions 
can thus be written 
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where c ruai ^,a ( c r a k &) ^ s tne destruction (creation) 
operator for a fermion with spin a, momentum k rel- 
ative to the Fermi point n ■ K{a{), and belonging to 
band ai. In ([!]) the sum over momenta is restricted to 
\k + q\, \k\, \k' — q\, \k'\ < A, and r; = ±1 indicates the 
branch of the Fermi surface where the particle resides. 
Conservation of crystal momentum is enforced by the 
function Sp. 

The allowed values of the couplings g are, of course, 
constrained by symmetries such as conservation of crystal 
momentum in the longitudinal (x) direction, spin conser- 
vation, parity under reflection about the x axis, (in the 
PBC case, also transverse momentum conservation), time 
reversal and particle-exchange symmetry These symme- 
tries severely restrict the allowed couplings, especially if 
there is no particular relation between the Fermi mo- 
menta at the given value of the particle density. 

In the OBC case, the allowed couplings can be clas- 
sified like for the two-legs system, [[lOj where two band 
indices a/ in (Eh belong to one band and the other two 
to another or to the same band, alternatively with paral- 
lel or antiparallel spin combinations, whereby the r; are 
fixed by x momentum conservation. ||llf On the contrary, 
for PBC the situation is more complicated since (a) some 
processes are forbidden by transverse momentum conser- 
vation, and (b) the degeneracy of the two bands a = 2 
and a — 3 allows for anomalous processes in which par- 
ticles are scattered among all three bands. In addition, 
more couplings are present at special particle densities 
both in the OBC (for example, at half filling) and in the 
PBC case. § 

The RG flow is generated by reducing the momentum 
cutoff A by infinitesimal steps and by modifying (renor- 
malizing) the parameters of the Hamiltonian in such a 



way that the partition function remains fixed. |9) At one- 
loop order, the differential renormalization of the param- 
eters is evaluated by integrating the momenta of the in- 
ternal lines over an infinitesimal shell of width AA at 
the cutoff in the corresponding diagram. § The differ- 
ential renormalization Ag a . a > (r±, ay, r$, 0:3 1 7-4, 0:4; r?,, a-i) 
of a given coupling is simply obtained by evaluating the 
"bubble" diagrams of Fig. 1 with zero external momenta 
and wave vectors, whereby the internal wave vectors q are 
restricted to the infinitesimal "shell" A — AA < \q\ < A. 
The integration is readily evaluated yielding from each 
diagram the contribution r)5 rs +r 6 g A g B g^r^j+f^j) T > 
where w 5 and v§ are the Fermi velocities of the inter- 
nal lines 5 and 6, respectively, and 77 is a symmetry and 
sign factor. The total contribution to the renormal- 
ization of a given couplings corresponding to given ex- 
ternal "legs" 1, 2, 3, 4 of the diagram is then obtained by 
summing over all possible sets of labels {ai = 1, ... ,3, 
ri = ±1, and a = ±1) for the internal lines 5 and 6 cor- 
responding to allowed couplings g A and gs for the points 
A and B. In particular, it follows as usually that cou- 
plings with n = r2 = r 3 = r^ do not renormalize at the 
lowest order. |]i],^|Jl0| The effect of the renormalization of 
Fermi velocities and momenta can be neglected at this 
order. 
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FIG. 1 . One-loop diagrams contributing to the renormal- 
ization of the couplings g(l;3|4; 2). 

After carrying out the burdensome but straightforward 
summation for the different couplings, by means of a 
computer algorithm for symbolic calculations, we end up 
with a set of quadratic RG differential equations, (40 
equations at half filling, 27 for general density, in the 
OBC case, 24 in the PBC case) that can be written: 



d 9i 
d t 



= A 
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where r = — In A, the indices and I label the allowed 
processes, and the matrix elements A*- 1 depend on the 
Fermi velocities. 

The precise expression for the system of differential 
equations (||) would be too long to reproduce here, 
however, it can be readily solved numerically by inserting 



initial values for the couplings depending on the particu- 
lar model. The purpose here is to determine which cou- 
plings are most relevant in the RG sense, that is, which 
ones of the couplings in (|l]) are important for the low- 
energy properties of the system. This is achieved by 
looking at the behavior of their RG flux. It turns out 
that, depending on the initial conditions, some of the 
couplings will diverge at a certain cutoff scale t = tq and 
some not (only in few cases, all couplings converge). In 
the weak-coupling limit the renormalized Hamiltonian, 
whereby only the relevant couplings are retained, can be 
constructed by fitting the asymptotic behavior of each 
coupling by the power law g t (r) « ( T p!l? T ) > an d by replac- 
ing each gi by its associated <?io in the Hamiltonian. H 
This procedure is controlled in the weak-coupling limit 
since one can always choose an arbitrary small U such 
that the asymptotic behavior occurs in the weak-coupling 
regime. 

In order to understand the low-energy physics that fol- 
lows from the renormalized Hamiltonian, we write it in 
terms of bosonic fields, |S| by expressing the fermionic 
destruction operators in the real-space representation as 



c r ,a,u(x) » C 'VX exp[i (r K (a) x - (r(f> p . a (x) - 6 p<a (x)) 
+o-(r(/)a, a (x) - 6<t,<x{x)))] , (3) 

where C is a constant of order unity, and 4> p , a (x) , 9 p , a (x) 
are dual canonical bosonic fields associated with the 
charge (p = p) or spin (p = a) modes on band a. 




FIG. 2 . Phase diagram for the three-legs Hubbard ladder 
with open transverse boundary conditions (OBC) as a func- 
tion of transverse hopping t p > 0, and particle density p > 1. 
Thick lines represent regions (phase lines) where a phase oc- 
curs for a strict relation between the Fermi momenta in the 
c7 — > + limit. Phase boundaries are indicated by dashed lines 
for the extended regions and by dots for the phase lines. Con- 
ventions for the labeling of the different phases are explained 
in the text. 

The resulting bosonic Hamiltonian consists of 
two kinds of couplings. Couplings of the density - 
density type (containing operators of the type 



c l, a A x ) c rM,a{x)cl, a , <a ,{y)cr' , a ',cr'(y)) contribute to the 
quadratic ungapped part of the Hamiltonian by mixing 
the different band modes, and renormalizing the Fermi 
velocities. The model restricted to density-density type 
couplings conserves spin and charge at each point of the 
Fermi surface and is the analogue of a Fermi liquid in 
quasi-one dimension, it is exactly solvable and shows 
spin-charge separation as well as a power-law behavior of 
correlation functions. However, the rest of the couplings 
(which will be here referred to as boson-interaction cou- 
plings) add to the Hamiltonian complicated "cosinus"- 
like terms that cannot be treated exactly. The standard 
way to deal with these terms, whenever they turn out to 
be relevant, is to expand the bosonic fields about one of 
the (equivalent) minima of the interacting Hamiltonian 
thus yielding a mass gap for the associated bosonic field. 

in 

We first analyze the region where all three bands cross 
the Fermi surface [region (A) of Fig. 2 ] and no spe- 
cial combinations of Fermi momenta occur. As it is clear 
from the discussion above, we are interested in deter- 
mining when some of the three charge and three spin 
modes (one for each band) acquire a gap due to the rele- 
vant processes. For this reason, we only need to consider 
the relevant couplings which are of the boson-interaction 
type. In this case, one has two "backward" couplings 
gb(i) reversing the spin in the two opposite branches of 
each of the two bands i = 1 and 3, and two processes 
5/b(l,3) transferring particles that move in opposite di- 
rections from band 1 to 3 and vice-versa (with different 
spin symmetry). Band 2 turns out to be asymptotically 
free. The associated boson-interaction term, obtained by 
using Eq. (^) reads: 

5b(l) cos(40 CTj i) +5b(3) cos(4ov, 3 ) 

+20 /6 (l,3) sin(20 ffil ) sin(2^, 3 ) cos[2(0 p ,i - 9 Pt3 )} . (4) 

Expansion about the minimum at 4>a,a = tt/4 (a = 
1,3) and 6 P - = 8 Ptl — 9 p , 3 = tt/2, yields a mass gap 
for the two spin modes </> CTi i and <j) a ,3, and for the odd 
charge mode 9 p -- Using the same notation as in Ref. 
||, whereby a phase with n gapless charge and m gap- 
less spin modes is referred to as CnSm, we have here a 
C2S1 phase (see Fig. 2 ). Note that, unlike the two- 
legs system, where a C1S0 phase occurs, || there are 
ungapped spin modes about half filling, which implies, 
for example, that the spin correlation function does not 
decay exponentially. 

The situation is quite different for the PBC system. In 
this case, several boson-interaction couplings can become 
relevant, including some of the processes that involve par- 
ticles in the three different bands. Two regions can be 
distinguished: When the Fermi velocity vfi of the first 
band is larger than the velocities of the other two (de- 
generate) bands «F23, the relevant couplings produce a 
C2S1 phase similarly to OBC. On the other hand, when 
vfi < VF23 only the uniform charge mode is ungapped 
yielding a C1S0 phase. Notice that this last region in- 
cludes half filling, where a spin-liquid phase is expected, 
due to the frustration of the periodic system. 

At half filling, a number of additional couplings are 
allowed for the OBC system. Due to the relation 



K(l) + K{2>) = 2K{2) = 7r, one can have umklapp pro- 
cesses with several combinations of Fermi momenta all 
equal to 2n. It follows that, in addition to the usual 
umklapp term for band a = 2 one has combined um- 
klapp processes for bands a = 1 and a — 3 as well as 
processes that scatter particles from and to the three dif- 
ferent bands. The behavior of the RG flow for these 
couplings turns out to be different depending on the 



value of t p . For t p /t > x c 
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Hamiltonian is restricted to the combined 1 — 3 umk- 
lapp terms (equal to <?/&( 1,3)) and to the interaction 
terms (|J) (and to unimportant density-density terms). 
Therefore, in addition to (^) one has a bosonic interac- 
tion of the form 5/5(1, 3) cos[2(0 Pi i + </)p i3 )][cos(2# Pi _) + 
2 cos(20 (7i x) cos(2</> CT ,3)] which leads to an additional gap 
in the p ,i + <j) P ,3 charge mode and to a C1S1 phase. Re- 
markably, the umklapp process on band 2 is suppressed 
by the interaction with bands 1 and 3 and is thus ir- 
relevant in this region. The system is thus a conductor 
(the current being carried by the band a = 2) with an 
ungapped spin mode. [fl3[ 

For tp/t < x c a number of additional couplings, in- 
cluding some of the processes that involve particles in 
the three different bands become relevant. The analy- 
sis of the complicated boson-interaction Hamiltonian Jll]] 
leads to a C0S0 phase with completely gapped charge and 
spin modes (like for the two-legs system at half filling 
||). One therefore has a metal-to insulator transition 
for decreasing t p at t p /t — x c . Notice that, regarding 
the spin modes, our results are consistent with the ex- 
periments (which probe an essentially isotropic system, 
thus t p /t w 1). H|l3|] It would be interesting to verify 
whether in anisotropic spin ladders (t p < t) a spin gap is 
present at half filling or whether our result is due to the 
weak-coupling condition. 

The half-filling condition for band 1 is fulfilled on the 
upper thick line of Fig. 2 . In this region a single gapped 
charge mode occurs. || We shall neglect other special 
cases which can arise due to particular combinations of 
Fermi momenta or when a band edge touches the Fermi 
level. § 

In region (B) of the phase diagram only two bands 
cross the Fermi surface. The analysis of this region is 
identical to that for the two-legs ladder || yielding the 
same set of RG equations and similar phases. Finally, 
in the remaining regions (C) and (D) only one band is 
involved and one recovers the results of the purely one- 
dimensional system. JsJ 

In order to check the dependence on the detailed fea- 
tures of the particular model, we have tested our results 
by including, alternatively, a nearest-neighbor U nn and a 
transverse U p interaction. While values of U nn /U up to 
0.5 essentially don't modify the phase diagram, it turns 
out that already with U p /U > 0.2 charge and spin exci- 
tations tend to be suppressed and additional gaps arise 
especially in the neighborhood of the phase boundaries, 
both at half filling and at finite doping. 

In conclusion, we have investigated the weak-coupling 



phase diagram of three coupled one-dimensional elec- 
tronic systems modeled by a Hubbard Hamiltonian (al- 
lowing eventually for some extensions) with open bound- 
ary conditions in the transverse direction. Consistently 
with experiments on S^CuaOs, spin modes are ungapped 
in a large region at finite doping and at half filling for 
tp/t greater than a critical value x c . In the same region, 
the system is a conductor for U small enough and a Mott 
transition occurs at half filling by decreasing t p . We have 
also analyzed the system with periodic transverse bound- 
ary conditions which shows, in contrast, a spin gap in a 
large region at and close to half filling. 

As it was clear for the two-legs case, [[| the interest of 
this analysis applied to a quasi-one-dimensional system 
relies upon the non-trivial phase diagram which arises 
already at weak coupling. Here, the system is strongly 
correlated even at infinitesimal U . By increasing U , one 
should expect a smooth deformation of the phase bound- 
aries and a widening of the isolated phase lines in Fig. 
2 , such that the general features of the phase diagram 
should not change dramatically for finite U (up to a cer- 
tain value |p^[| ). 

We thank G. C. Strinati, B. Brendel, J. Voit, H. J. 
Schulz, and W. Hanke for stimulating discussions and 
precious suggestions. The author acknowledges research 
support from the ECC under the network program N. 
ERBCHRX-CT940438. 



[1] M. Fabrizio, A. Parola, and E. Tosatti, Phys. Rev. B 46, 
3159 (1992); M. Fabrizio, Phys. Rev. B 48, 15838 (1993). 

[2] S.R. White, R.M. Noack, and D.J. Scalapino, Phys. Rev. 
Lett. 73, 886 (1994). D. J. Scalapino, Nature 377, 12; 

[3] L. Balents and M. P. A. Fisher, Phys. Re v. B in press. 

[4] H. J. Schulz, preprint |co~nd-mat 941209^ (1994). 

[5] S. Gopalan, T. M. Rice, and M. Sigrist, Phys. Rev. B 49, 
8901 (1994). 

[6] K. Kojima, et al, Phys. Rev. Lett. 74, 2812 (1995). 

[7] V. J. Emery, in Highly- Conducting One-Dtmensional 

Solids, edited by J. T. Devreese, R. E. Evrard, and V. E. 

van Doren, Plenum Press, New York (1979). 
[8] See, e. g., J. Solyom, Adv. Phys. 28, 201 (1979). 
[9] See, e.g., K. G. Wilson and J. Kogut, Physics Reports 

12, 75 (1974). 

[10] C. M. Varma and A. Zawadowski, Phys. Rev. B 32, 7399 
(1985). 

[11] Details will be given elsewhere: E. Arrigoni (unpub- 
lished) . 

[12] Despite the presence of the boson-interaction term one 
can still obtain an ungapped mode whenever a dual sym- 
metry cf> <-> 8 is present (see Ref. this situation, how- 
ever, never occurs in our results. 

[13] One expects that the metallic half-filling phase (for 
t/tp > x c ) becomes insulating at a certain finite U — U c 
due to the relevance of the umklapp scattering in band 2 
at higher order. On the contrary, spin modes should be 
enhanced by increasing U. 



